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We analyze the stabilization of unstable steady states by delayed feedback control with a periodic 
time-varying delay in the regime of a high-frequency modulation of the delay. The average effect 
of the delayed feedback term in the control force is equivalent to a distributed delay in the interval 
of the modulation, and the obtained distribution depends on the type of the modulation. In our 
analysis we use a simple generic normal form of an unstable focus, and investigate the effects of 
phase-dependent coupling and the influence of the control loop latency on the controllability. The 
main advantage of the variable delay is the considerably larger domain of stabilization in parameter 
space. 
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I. INTRODUCTION 

The possibility to stabilize unstable periodic or sta- 
tionary states embedded in chaotic attractors has been 
elaborated more than two decades ago [l| . The main idea 
consists of using small external perturbations to force 
the system to follow one of its stable manifolds. These 
perturbations are applied at specific instances when the 
chaotic trajectory is close to the desired periodic orbit. 
This seemingly straightforward theoretical concept has 
caused a revolution in applied nonlinear science. It has 
been realized that the control of chaos could have a sig- 
nificant outcome in real-world experiments, where one 
could generate different kinds of ordered behavior from 
an utterly erratic one [^ ■ 

A conceptually simple method to stabilize unstable 
equilibria and periodic orbits is the time-delayed feed- 
back control (TDFC) introduced by Pyragas 0. Here, 
the perturbation has the form of a continuous feedback 
constructed from the difference between some suitable 
scalar signal obtained from the system and the same sig- 
nal delayed by a constant time r. The difference signal is 
amplified by the gain factor K and then re-injected into 
the original system. For a certain choice of the feedback 
gain K and the delay time r, the control of the unstable 
state can be realized, in which case the feedback force 
vanishes by construction, making the method essentially 
noninvasive. 

Since a detailed knowledge of the target state is not 
required and the controller is very robust with respect 
to noise, the Pyragas control has been widely used in 
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control experiments [j|, and substantial work has been 
done to understand the control mechanism analytically 
[5l4lOl| . The method was subsequently improved by intro- 
ducing multiple delayed signals into the feedback control 
force, both with commensurate or incommensurate de- 
lay times [lll - ll4l | . Another improvement of the method 
was achieved by introducing a time-varying delay into 
the delayed argument of the control force which could 
be realized experimentally by changing some character- 
istics of the delay line in a deterministic (periodic) or 
a stochastic fashion [l5l - [l7t . For example, in coupled 
laser systems the delay time could be modulated by pe- 
riodically changing the distance of the lasers or external 
resonator, e.g. by piezoelectric modulation. For systems 
acting on a slower timescale, e.g. electronic or mechanic 
devices, one may take advantage of digital delay lines in- 
stead, where the time-functional form of the delay modu- 
lation could be controlled by an external clock frequency 
modulator [l8|. In this way, one could practically re- 
alize the variable- delay feedback controller and stabilize 
unstable steady states and periodic orbits over a much 
larger domain of control parameters. To keep the method 
noninvasive also in case of unstable periodic orbits, the 
delay modulation should be chosen appropriately, e.g., 
in a form of a (non)pcriodic change of the delay between 
multiples of the basic period of the uncontrolled unstable 
orbit [13. 

In this paper, we aim to provide a deeper analyti- 
cal insight into the mechanism of stabilization of un- 
stable equilibria by the delayed feedback method with 
a time-dependent delay time. A general description of 
the variable-delay feedback controller is given in Sec. II, 
where it is shown that the method could be explored an- 
alytically in the domain of a high-frequency modulation 
of the delay. In this case, we use the formalism of dis- 
tributed delays [20, [21| , where the average contribution 



of the time-varying delay is represented by an integral 
kernel describing a particular delay distribution. In Sec. 
Ill, we consider a phase-dependent coupling of the control 
force, which is a particular non-diagonal generalization of 
the standard diagonal coupling scheme. As a model sub- 
ject to control we use a generic two-dimensional linear 
system that has an unstable steady state of focus type. 
This model mimics the two-dimensional center manifold 
of a general non-linear system, capturing the essence of 
the dynamics in the vicinity of its unstable fixed point. 
In Sec. IV, we investigate the influence of control-loop 
latency on the efficiency of the controller. We conclude 
our findings in Sec. V. 



compared to the uncontrolled system's dynamics de- 
scribed by an intrinsic frequency w, the linearized system 
([2]) with a modulated delay has the same asymptotic sta- 
bility properties as the averaged linear distributed-delay 
system 

x(t) = Ax(t) + K ( / pi9)x{t -9)d9- x(i) ] . (4) 



The probability density function p{0), i.e. the 
distributed-delay kernel, is defined in a way that p{9)d9 
gives the fraction of time for which r(t) lies between 9 
and 9 + d9, satisfying p{9) > and the probability nor- 
malization condition 



II. VARIABLE-DELAY FEEDBACK CONTROL 

We consider a general n-dimensional dynamical system 
X = r(x), where F is the vector field describing the sys- 
tem's dynamics, and x = x(f) e M" is the state column 
vector. The equilibria x* of the system are the solutions 
of F(x* ) = 0, and the stability of a particular equilibrium 
X* is determined by the eigenvalues of the Jacobian ma- 
trix A = 9xF(x) calculated at x = x*. In the following 
we assume that x* is an unstable equilibrium, having at 
least one eigenvalue with a positive real part. 

Under a linear variable-delay feedback control, the 
original system is transformed into 



i(t) = F(x(t))+K[x(t-r(i))-x(t)] 



(1) 



where K is an n x n feedback gain matrix, and T(t) is 
the time-dependent delay time. By choosing K and T(t) 
appropriately, we aim to stabilize the unstable equilib- 
rium X* . Without loss of generality, we assume that the 
equilibrium has already been moved to the origin by a 
translation of axes, such that x* = 0. We thus linearize 
the controlled system in the neighborhood of the origin 
to obtain 



x(t) = Ax(t) + K [x(i - T(t)) - x(t)] 



(2) 



We consider a deterministic modulation of the vari- 
able time-delay T{t) in a fixed interval around a nominal 
(average) delay value tq 



T{t)=T^+ef{vt), 



(3) 



where / : R — >■ [— 1, 1] is a 27r-pcriodic function with zero 
mean, and e and u are additional variable-delay control 
parameters denoting the amplitude and the (angular) fre- 
quency of the modulation, respectively. The stability of 
the origin can be inferred by numerically integrating the 
linear variable-delay system ([2]) for different values of the 
control parameters K, tq, e and v, thus determining the 
domains in the parameter space for which the stabiliza- 
tion can be realized. 

In the limiting case of a high-frequency modulation 
[la . |22| | , when the parameter v becomes sufficiently large 



p{9) d9 = 1. 



(5) 



In practice, the transition to the distributed-delay limit 
case does not require very large modulation frequencies, 
and therefore variations of the delay in the experiment 
are a very practical way to create different types of delay 
distributions. 

For a continuous modulation in an e interval around 
Tq in form of a sawtooth- wave (triangular) modulation of 
the delay. 



T{t) = To + e 



^ , vt 

2 — mod I 

27r 



1 



(6) 



p{9) is a constant in the interval of the delay variation. 
Under the probability normalization condition ([5]) we ob- 
tain a uniform distribution 



p{o) = l; 



, 9 e [To-e,To + e], 
, elsewhere, 



(7) 



which does not depend on the skcwness of the sawtooth 
function. For a sine-wave modulation, 



r(<) = To + £ sin(t/f). 



(8) 



the time interval dt during which the delay r changes by 
dr is given by 



dt^ 



dr 



dr 



i^ecos 



ii^t) y^e'' - (r - ro)2 



(9) 



The fraction of time dt within a half-period ix jv of the 
delay function t is equivalent to the product p(T)dT. In 
terms of our previous notation, we get 



P(^) 



TTv/e^ ~{9- ro)2 



(10) 



A periodic change of T(t) between two fixed values tq — £ 
and To 4- £ for the same time duration. 



r(t) = To +£sgn[sin(i/i)]. 



(11) 



results in a square wave (rectangular) modulation with a 
two-peak distribution 



p{e) = \ (6(9 - To + e) + Sie -To- e)) 



(12) 



In an analogous way, we may obtain the corresponding 
distributed delay kernels for other types of delay modu- 
lations. 

Applying the exponential ansatz x(i) ~ exp(At) in Eq. 
(HI, we obtain the characteristic equation for the eigen- 
values A 



det 



AI-A 



1 



pi0)e 



^Ae, 



de] K 



0, (13) 



where I is the identity matrix. Since p{9) is nonzero only 
between tq — e and tq -I- £, we have 



det 



AI A 



1-e 



'Ato 



x(A,e))K 



= 0. 



The quantity 



x(A,e) 



+e 



p(to + 9)e-'''d9 



(14) 



(15) 



summarizes the effect of a given modulation, and for the 
above modulation types reads 

{ ^™.^ '^•^ , sawtooth- wave, 
Io{Ae), sine-wave, (16) 

cosh(A£), square- wave, 

where /q is the modified Bcssel function of the first kind of 
order 0. In the non- modulated case, x(A, 0) = 1, we have 
the usual TDFC. For general distributed delay, the delay 
term exp(— Atq) is replaced in the characteristic equation 
(IT51) by the Laplace transform of the delay kernel p{9) 

The function x(A, e) also allows for a compact explana- 
tion of the mechanism of variable-delay feedback control 
(VDFC). The expressions listed in Eq. p^ have regions 
{Ae e C} for which |x| < 1, mostly for Re(A) < Im(A). 
In particular, values of x « can occur, meaning that the 
delay term effectively has been suppressed and vanishes 
from the characteristic Eq. (J14p . In the original equations 
of motion Q this situation can be interpreted as de- 
structive interference of the delay signals covered by the 
distribution. Positive and negative phases of the spiral 
oscillation in the neighborhood of the steady state cancel 
each other out creating a reference term x^oi{t) ~ x* for 
a simplified controller 



i(t) = Ax(i) + K (x,ef(i) - x(i)) 
« Ax(i)+K(x*-x(i)) . 



(17) 



The same idea can be regarded as the motivation for the 
original TDFC method. The advantage of VDFC lies 
in an improved approximation of the target state by the 
delay terms. Fig. [T] illustrates the control mechanism 




FIG. 1: (Color online) Sample of a trajectory x{t) (black) 
in the vicinity of a steady state (red/gray horizontal line) of 
a system subjected to variable-delay feedback control with a 
fast sawtooth wave. The corresponding uniform delay distri- 
bution ([7} creates a sliding average between t — to — e and 
t — To + e (area highlighted in blue/gray) which at time t re- 
sults in the reference signal (blue/gray, oscillating). The con- 
trol force is constructed from the difference signal as marked 
by the vertical black arrow. 



for an almost ideal situation, in which the steady state 
is approximated by the averaged delay signal. The in- 
stantaneous part of the coupling term stabilizes the sys- 
tem while pulling it towards the reference signal x,-ci{t) 
which due to the small x-value has a smaller amplitude 
than the remaining instantaneous x(t). The resulting 
stabilization of the steady state shows a high robustness 
against parameter detuning, because the simplified con- 
trol mechanism is insensitive to the phase relation be- 
tween actual signal and reference signal and also works 
for a wide range in the coupling gain K. However, the 
described scenario relies on a small x- value, which in gen- 
eral is not trivial to obtain. The delay terms can almost 
never be suppressed completely, so that only a rigorous 
consideration of the characteristic equation (|14p reveals 
the full capability of the control method. Eq. ([T4|) is tran- 
scendental with respect to A, possessing an infinite set of 
complex solutions {Ai} £ C. The origin can be stabilized 
for those values of the control parameters (K,to,£) for 
which all the eigenvalues {Ai} have negative real parts. 
The stability domain in the parameter space (K,ro,e) 
can be calculated numerically given the modulation dis- 
tribution. Note that the control parameter v is lost in 
the transition to distributed delays. 

As in the usual TDFC and extended TDFC control 
schemes, the presence of torsion is necessary for the pro- 
posed control method to be able to stabilize equilibria. 
Torsion means that in its unstable subspace the steady 
state is only of the spiral type. To show this property, 
we consider the characteristic quasipolynomial -ff(A) = 

det 



AI - A + (1 - e-^^«x(A, £)) K 



It is easily shown 



that this quasipolynomial is positive for A -^ oo, and for 



A = reads H{0) = det -A = 



nLi(- 



s„), where e. 



arc the eigenvalues of A. If A possesses an odd number 
of positive real eigenvalues, then H{0) < 0, and there ex- 
ists at least one positive real root of H{A) = 0, meaning 
that the fixed point cannot be stabilized by the proposed 
method. Although the analysis is done in the distributed- 
delay limit, numerical simulations show that this odd- 
number limitation persists also for low-frequency modu- 
lations of the delay r(t). 

The results of this section will be exploited in the 
following to perform a comparative analysis between 
the variable-delay feedback control in a distributed-delay 
limit and the standard delayed feedback control by us- 
ing a simple normal form model as a representative of 
a large class of nonlinear dynamical systems. Specifi- 
cally, we will investigate the effects of including a phase- 
dependent coupling matrix and the influence of nonzero 
latency times in the control scheme. 



III. PHASE-DEPENDENT COUPLING 

We will apply the variable-delay feedback control to an 
unstable steady state of focus type. This system repre- 
sents a generic model of an unstable steady state slightly 
above a Hopf bifurcation. In center manifold coordinates. 



J 



the dynamics of the system is given by Eq. ([2|), where 
X = {x, y)^ is the two-dimensional column vector, and A 
is a 2 X 2 matrix 



A cj 

— Ld A 



(18) 



that describes the dynamics in the absence of control. 
Since the stability of the free-running system is deter- 
mined by the eigenvalues A ± i oj of A, choosing A > 
and a; 7^ we model an unstable focus at the origin. 

In this section, we investigate a phase-dependent cou- 
pling of the control force, which is relevant, e.g., in sta- 
bilization of laser devices where the optical phase occurs 
as an additional degree of freedom [23, [2J1. Such cou- 
plings have also been used in the context of refuting the 
odd-number limitation of delayed feedback control in au- 
tonomous systems [^, llflj l2a]; ^o anticipate chaos syn- 
chronization [2g], and more recently, to control different 
synchronous oscillatory states of delay-coupled oscillator 
networks [231 ■ The phase-dependent coupling is realized 
via a rotational matrix K containing a variable phase (/?, 
and it enters as an additional multiplicative factor to the 
control force 



x{t) 

m 



A LU 

—uj A 



x{t) 



cosip —smip\ f x{t — T{t)) — x{t) 



suvLp costp J \y{t~T{t))-y{t) 
where K is the feedback gain. In the distributed-delay limit, the system becomes 



x{t) 

m 



A LO 
— LO A 



x{t) 

yit) 



, , , cosi^ — siniysX / Jq p{9)x{t — 9) dO — x{t) 
^ ' sin ^ cos ^ j 1^ /;° p{e)y{t - 0) dO - y[t) 



and the stability of the origin is determined by the roots A of the characteristic equation 

[A - A + A' (1 - e-^^'^'xiK e)) cos Lpf +[uj + K{l- e-^"'«x(A, e)) sin y^] ^ = 0. 

I 



(19) 



(20) 



(21) 



This equation can be further simplified to 

A + Ae^"^ (1 - e~^^«x(A, e)) = A ± iu. (22) 

The control method is successful if there exists a 
nonempty set in the four-dimensional parameter space 
(A, To, e, (p) for which the real parts of all the character- 
istic roots A are negative. 

The domain of control in the plane parametrized by 
the feedback phase ip and the feedback gain K is given 
in Figs. [5] and [31 The grayscale (color code) indicates 
only those control parameters (ip, A) for which the lead- 
ing eigenvalues have negative real parts, and the control 
is more robust as these values are more negative. Panels 
(a)-(d) in Fig. [5] depict control domains for a sawtooth- 
wave modulation corresponding to a fixed nominal delay 



To = 1 and different modulation amplitudes e: (a) e = 0, 
(b) e ~ 0.3, (c) £ ~ 0.5 and (d) e — 0.9. The parameters 
of the unstable focus are fixed at A = 0.1 and u = it. 
The chosen nominal delay tq = 1 is an optimal value 
in the non-modulated control case and tp = Q (diagonal 
coupling). From the corresponding control domains, it is 
observed that an increase of the modulation amplitude 
£ leads to an enlargement of the control domain in the 
direction of the positive (/3-axis, thus destroying the sym- 
metry of the stability island centered at (p = 0. A similar 
result is observed for a sine-wave modulation in panels 
(a)-(b) in Fig. [31 corresponding to e = 0.5 (a) and e = 0.9 
(b). A different behavior is observed in panels (c)-(d) in 
Fig. [31 depicting the corresponding stability domains for 
a square-wave modulation. In the latter case, the stabil- 




FIG. 2: (Color online) Control domain in the [ip, K) plane for 
a sawtooth-wave modulation of the delay for different mod- 
ulation amplitudes: (a) e = 0, (b) e = 0.3, (c) e — 0.5 and 
(d) e — 0.9. The nominal delay is fixed at tq — 1, which is 
optimal value in the non- modulated case (TDFC). The pa- 
rameters of the free-running system are A = 0.1 and uj = n. 
The grayscale (color code) depicts only those values of the 
control parameters for which the largest real part of the com- 
plex eigenvalues A is negative, indicating a successful control. 
The depicted solid lines that enter into a description of the 
stability boundary are calculated from Eqs. (|26|) and (|27|l . 
and the dashed lines from Eqs. (|26p and (|28[1. For clarity of 
the picture, we depict only a few branches of the boundary 
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FIG. 3: (Color online) Same as Fig. [2] for a sine-wave mod- 
ulation of the delay: (a) e = 0.5, (b) e — 0.9, and for a 
square-wave modulation: (c) e — 0.5, (d) e — 0.9. Other 
parameters are the same as in Fig. [S] 



ity domain increases more rapidly as e is increased from 
zero, until e reaches a critical value after which the do- 
main decreases, resulting in an instability stripe at ip ~ 
when e = 1 . This non-monotonic change of the stability 
domain is typical for a square-wave modulation, and it 
has been observed in other circumstances. 

In order to understand the mechanisms behind these 
numerical results, we first recall the intuitive argument 
about destructive interference and the simplified con- 
troller from the previous section. From Eq. ((20)) we di- 
rectly see that the stabilizing diagonal elements of the in- 
stantaneous control term are weighted with cos f. Since 
the onset of stability /Cmin is dominated by the com- 
petition between A > and —Kcosip < 0, we ex- 
pect Kynin oc X/ COS if, scc also the detailed reasoning in 
Eq. ([50)1 below. For \(p\ > 7r/2 the instantaneous part of 
the coupling becomes repelling. Therefore with regard to 
the depicted control mechanism we expect the control to 
largely fail in this region. All simulations clearly support 
this picture. Small exceptional regions of stability ex- 
ceeding this boundary arise from the non- vanishing delay 
terms which can support stabilization given an optimal 
parameter constellation. However, for large values of the 
mean delay tq this type of support fails and \ip\ = 7r/2 
becomes a strict limit of control. 

The asymmetry can be explained by the interaction of 
the non-diagonal elements of system terms and control 
terms in Eq. (pn|) . If cj were zero, the control term would 
have no preference for positive or negative values of f 
because the change of tp -^ —ip would be equivalent to 
mirroring the complete system at one coordinate axis, 
e.g. X ^ —X. A non-zero value of the spiral frequency 
to breaks the symmetry and leads to different resulting 
frequencies fl = Im(A) under control for different signs 
of if. From the interference mechanism one can conclude, 
that at least for the continuous distributions ([7]) and pU]) 
high frequencies are more favorable for successful control, 
because a low x-value corresponds to many oscillation 
periods covered by the distribution. Indeed, if we have 
a closer look at the parameter constellation underlying 
Fig. [21 we see that for negative (p the imaginary part of 
the most unstable mode tends to vanish, thus inhibit- 
ing the interference effect and leading to control failure. 
In contrast, for positive ip the frequency increases with 
increasing K , thus improving control efficiency. A para- 
metric plot of different A^ in the complex plane reveals 
this property, see Fig. [H 

A detailed analytical investigation of the control do- 
main can be done by using the characteristic equation 
(|22p and noticing that the transition from instability to 
stability occurs at the boundaries of the control domain, 
where the leading eigenvalues A are purely imaginary. 
Setting A ~ in in Eq. ((22)) . and separating real and 
imaginary parts results in a system of two real-valued 
equations 

A' cos p — Kxiin, e) cos(Oro ± </?) ~ A, (23) 
fl^ Ksinp: + Ksin{i}To±(p)x{'i'^,£) = i"^, (24) 
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FIG. 4: (Color online) Characteristic exponents A in the 
complex plane parametrized by the feedback gain K. Left: 
ip = —1, right: (p — +1. Depicted are five leading eigenvalue 
curves for K £ [0, 6]. The gray-filled circle (yellow online) cor- 
responds to K = and the black dots correspond to K = 6. 
The eigenmodes created by the variable-delay control origi- 
nate from Re(A) = — oo. Other parameters as in Fig. [2};. 



which is an implicit parametric representation of the 
boundaries of the control domain in the parameter space 
{K,To,e,ip), in which the eigenfrequency V, has the role 
of a parametrization variable. We have taken into ac- 



J 



count that for the above mentioned distribution kernels, 
p{6) is even around the nominal delay tq, meaning that 
x{i^,s) is a real- valued function 



sin(n£ 



x{in,e) 



sawtooth- wave, 
sine- wave, 



cos (J7e) , square- wave 



(25) 



where Jq is the Bessel function of the first kind of or- 
der 0. Considering the positivity of the nominal delay tq 
and the feedback gain K, as well as the positivity of the 
parameters (/?, A and a; used in the current analysis, elim- 
ination of the phase parameter ip from the last system of 
transcendental equations leads to an expression for K in 
terms of the eigenfrequency D, 



K{n) 



A2 + (w - ny 



\nT„) + [x{in,e)-cos{nT„)Y 



(26) 



Taking into account the multivalued properties of the arc- 
sine function, one obtains in a similar manner the analyt- 
ical expressions for the phase parameter if in dependence 

on rt 



M^) - arcsin (^^(^)X(^»,e)[Asin(Oro)^+|c. -^»)^cos(»ro)] - {. - n)^ ^^^ 

M^) = - arcsin ( Kjn) ^^^"' '^ ^' ""^"""^ + ^T "^ T^"""^^ ' ^^ ' "^ )+ (2n + 1)., 
\ A-^ + (cj - 52)"^ / 

I 



(27) 
(28) 



where n is a non-negative integer. Equations 
describe the boundary of the control domain in Figs. 
[5] and [31 where the solid and dashed lines represent 
branches oi ipi{n) and (p2(f^), respectively. 



In the case of zero feedback phase ((^ = 0), the phase- 
dependent feedback force is reduced to a diagonal cou- 
pling. In the absence of modulation (TDFC), the domain 
of control in the plane spanned by the feedback gain and 
the nominal delay consists of several distinct stability is- 
lands centered at odd values of tq, separated by regions 
corresponding to even tq for which the stabilization fails 
for any K. This observation can be inferred from Eqs. 
dMI-dMl) setting 1^ = and £ = 0. When flro = {2n+l)TT, 
then K„iin = A/2 and il = to. On the other hand, when 
flTQ = 2n7r, the control fails for any finite value of K. 
A detailed analysis of the control boundaries in this case 
(TDFC) is provided in Ref. p. By the same reasoning, 
in the presence of modulation (e > 0) and zero phase. 



the corresponding minimum feedback gain values are 

T^— — ^. vT' To = (2r7. + l)7r/w, 

i^min(e) - { (29) 

■pj TT TT, To = 2m: /lu, 

leading to a reconfiguration of the control domain de- 
pending on the type of the delay modulation [l5| . Specif- 
ically, it is observed that the instability stripes at even 
values of tq cease to exist, and the stability islands are 
starting to interconnect as soon as e > 0. Although for a 
sawtooth and sine-wave modulation the enlargement of 
the control domain with increasing e is monotonic, an 
oscillatory behavior is observed for a square- wave modu- 
lation due to the form of the x function (i.e. cos(a;e)) in 
the denominator of the expressions for K^i^. In the lat- 
ter case, the stability islands centered at tq = {2n+l)TT / uj 
are first enlarged and eventually interconnected with in- 
creasing £, up to some s value after which the control 
domain shrinks, gradually collapsing into several stabil- 
ity islands centered at tq equal to even multiples of tt/lo. 
For (fi j^ and a variable-delay control force (e > 0), 



the minimum feedback gain at the above speeific values of 
To cannot be deduced from Eqs. ([221)- dMl) by following 
the previous lines of deduction. Nevertheless, we can 
make approximate predictions for small values of (p, when 
r^To is an integer multiple of tt, and we have Vl k, lo from 
Eq. (|24l). Hence, in the regime of small values of (p, the 
minimum feedback gain is 



A 



Krain[e.V) 



[1 + x(iw,e)J cosip' 



X 



To w (2n + 1)tt/u, 



To ~ Iwn juj. 



[1 -x(«w,e)Jcos(^' 

The resulting expressions for iiTmin show that in the 
absence of any delay modulation and for small (ys, the 
minimum feedback gain at the "optimal" delay values 
To = (2n + l)7r/w is if mm = A/(2cos(/3), containing the 
dependence on (p via cos p in the denominator. Conse- 
quently, the principal stability island in ((^, IC) paramet- 
ric plane is centered at (y9 = 0, for which A'min = A (see 
panel (a) in Fig. [5]). On the other hand, at tq ~ ^n-K juj^ 
the stabilization cannot be achieved for any K^ regardless 
of the value of the phase parameter <~p. For a non-zero 
modulation amplitude e at small ip^ the values of iiTmin de- 
pend on the type and the amplitude of the delay modula- 
tion. Although the form of x(iw, e) for sawtooth and sine- 
wave modulations is such that the control in the variable- 
delay case is generally possible for all nominal delays 
''0 ^ £7 this is not the case for a square- wave modulation. 
Namely, in the latter case, ^(^'^iS) = cos(aje), giving a 
non-monotonic behavior of -ffmin with increasing modu- 
lation amplitude e. Specifically, for e = (2?! + l)7r/cj, the 
control fails at To = (2n + 1)t:/uj for any K , but it is 
optimal at To = 2m: /lo. 

To further investigate the effects of the modulation am- 
plitude e on the control efficiency, we have calculated the 
stability domains in the parametric plane {K, e) for dif- 
ferent modulation types and different values of the phase 
parameter ip. The results are depicted in Figs. [SJ H] and 
[7] for sawtooth-, sine- and square- wave modulation, re- 
spectively. Different panels in each figure correspond to 
different phases: (a) 1^9 = 0, (b) </j = 7r/8, (c) p) = tt/A, 
(d) p = 3n/8. The parameters of the unperturbed sys- 
tem are A = 0.1 and cj = tt, and the nominal delay is 
fixed at To = 1 as before. As the modulation amplitude 
increases from zero, a considerable reconfiguration of the 
stability domain is observed, depending on the type of 
the delay modulation. Increasing the modulation am- 
plitude enlarges the interval of K for successful control, 
although the enlargement is not necessarily monotonic, 
and the domain may consist of several disconnected in- 
tervals. The improvement of the delayed feedback con- 
troller from including time- varying delay is evident from 
the diagrams. For example, for a phase parameter value 
p = 37r/8 (panel (d) in Figs. [5l [6] and [T)) , while the orig- 
inal Pyragas method fails for any K, the variable-delay 
feedback method is able to stabilize the fixed point for 
certain values of the modulation amplitude e > 0. 




FIG. 5: (Color online) Control domain in the {K, e) plane for 
a sawtooth- wave modulation of the delay for different values 
of the feedback phase p: (a) <^ = 0, (b) p = vr/S, (c) p — 7r/4 
and (d) p = Stt/S. The nominal delay is fixed at tq = 1. The 
parameters of the free-running system are as in Fig. [2] 




FIG. 6: (Color online) The corresponding control domain in 
the {K, e) plane for a sine-wave modulation of the delay for 
different values of the feedback phase p: (a) (/? = 0, (b) p = 
7r/8, (c) p = 7r/4 and (d) p = 37r/8. The other parameters 
are as in Fig. [S] 



The parametric representation of the stability bound- 
aries in {K, e) plane can be obtained from Eqs. ((23)) and 



X(*^,e) 



A sin(fiTo + p) + {cj — il) cos(r2To + p) 
sin(OTo) 



(31) 



A sin p + {lo — Q) cos p 



A sin(ilTo + p) + {uj — Q) cos(f7To + p) 



(32) 




FIG. 7: (Color online) The control domain in the {K, e) plane 
for a square- wave modulation of the delay for different values 
of the feedback phase ip: (a) 99 = 0, (b) ip = vr/S, (c) tp = 7r/4 
and (d) ip = Svr/S. The other parameters are as in Fig. [S] 



In the second equation, e is implicitly contained in 
x(iJ7,e), which is a real function in the delay modula- 
tion cases considered above (see Eq. E5|) . In general it 



is not possible to invert the function x analytically. One 
should treat the last system as implicit parametric rep- 
resentation of the control boundaries and seek for the 
stability curves numerically. The boundaries of the cor- 
responding control domains are given in Figs. EHZlby the 
solid line. 



IV. CONTROL-LOOP LATENCY 

In an experimental realization of the control method, 
one has to take into account the latency of the feedback 
circuit due to the time required for the generation of the 
feedback control signal and its reinjection into the sys- 
tem. In laser systems, the latency is associated with the 
time the light takes to traverse the distance between the 
laser and the Fabry-Perot controller. It has been shown 
that latency time always acts destructively upon the con- 
trol domainSjreducing the effectiveness of the controller 

In this section, we investigate the effects of latency 
on the variable-delay feedback control in the distributed- 
delay limit. The latency time 5 enters as an additional 
constant time-delay in the control force, and the system 
now reads 



J 



i(i) 

m 



A LO 
—LU A 



x{t) 

y{t) 



In the distributed-delay limit, the system becomes 



m \ 
m 



A LU 
—LU A 



x{t) 

vit) 



K 



x{t -S- T{t)) - x{t - 5) 
y{t - 5 - T{t)) - y[t - 5) 






de - x{t - 6) 

de - y{t - 6) 



r 



(33) 



(34) 



leading to the characteristic equation 

[A - A + Ke-^^ (1 - e-^"«x(A, e))] ' + w^ = 0, (35) 
that can be further simplified into 

A + Ke'^^ (1 - e-^^«x(A, e)) = A ± Ilu. 



(36) 



Upon separating the real and imaginary parts, one ob- 
tains a system of two real-valued equations. We again 
consider the control boundary given hy A = iil which is 
given by the implicit parametric representation 

Kcos{n6)~Kcos{Q{S + To))xii^,£) = A, (37) 
n-Ksm{nS) + Ksm{n{S + To))x{in,e) = ±lu{38) 

Figure |8] depicts the domain of control in the {to,K) 
parametric plane for a sawtooth-wave modulation at a 
fixed modulation amplitude e = 0.5 and increasing la- 
tency time S: (a) S = 0.1, (b) S = 0.2, (c) S = 0.3, (d) 



6 — 0.4. The control domains for a sine-wave modula- 
tion corresponding to S ~ 0.2 and 5 = 0.4 are depicted 
in panels (a) and (b) in Fig. [HI and panels (c) and (d) 
are related domains for a square-wave modulation. The 
parameter values of the uncontrolled system are set to 
A = 0.1 and oj = tt as before. For a zero latency time, 
the control domain fills the whole depicted range of the 
parametric plane down to some minimum value of K, 
which for tq equal to an integer multiple of tt/lu is ex- 
plicitly given by Eq. ((29)) . It is observed that increasing 
latency time reduces the area and the robustness of the 
control domain. In analogy to Sec. II, one can derive 
approximate expressions for i^min for small values of the 




latency S 



To 



TO 



FIG. 8: (Color online) Control domain in the {to, K) plane for 
a sawtooth-wave modulation of the delay for different values 
of the latency time: (a) 5 — 0.1, (b) 5 — 0.2, (c) 5 — 0.3 and 
(d) 5 = 0.4. The modulation amplitude is fixed at e = 0.5, 
for which the stability domain fills almost the entire depicted 
part of the parametric plane in the case of zero latency. The 
depicted solid lines that enter into a description of the stabil- 
ity boundary are calculated from Eqs. (|40fl and (|4ip . and the 
dashed lines from Eqs. (|40p and (|42p . The parameters of the 
free-running system are as in Fig. [2] and (fi = 0. 



Kmini£,S) 



A 

[1 + x(*^7 £)J cos(cj(5) ' 

A 
[1 -x(iw,£)Jcos(w(5)' 



To = {2?i + l)n/oj, 



To — 2nTr/u!. 



(39) 

In parallel to the conclusions of the previous section, in- 
clusion of a variable time-delay in the feedback control 
force generally enlarges the control domain, making the 
control possible even for those values of tq for which the 
control always fails at any K in the non-modulated case. 



To describe the boundaries of the control domain in the 
(to,K) plane analytically, we algebraically manipulate the 
system ([571) - ([55)l and find two families of branches of so- 
lutions for K and tq parametrized by the eigenfrequency 



J 



K{n) 

T0,l(^) 

To,2m 



A cos{nS) ~{u}~n) am{nS) ± y/[\cos{nS) ~{u}-n) sm{nS)]^ - (1 - x^(*^, e))[A2 + (w - f7)2] 



1 
n 
1 



l-X^i^.s) 



arcsm 



— arcsm 



A sin(fi(5) + {lu -Q) cos(fi(5) 
K{n)x{tn,e) 
Asin(OJ) + (cj - ri) cos{il6) 



K{n)x{in,e) 



2mT 



+ (2n + l)7r 



(40) 
(41) 
(42) 



r 



where n is a non-negative integer that takes care of the 
different branches due to the multivalued arcsine func- 
tion involved in the boundary description. The control 
domain boundaries in Figs. [5] and [S] are represented para- 
metrically by Eqs. (|40)) - (|42|) . where the branches related 
to tq.i are given by the solid lines, and the branches of 
To, 2 with dashed lines. 

The influence of the modulation amplitude e on the 
stability of the fixed point can be deduced by observ- 
ing the control domains in the {K, e) parameter plane. 
Figures [101 HD and [12] depict the control domains for 
sawtooth-, sine- and square- wave modulation, respec- 
tively, for a fixed nominal delay tq = 1 and increasing val- 
ues of the latency parameter 8: (a) 5 = 0.1, (b) 5 = 0.2, 
(c) 5 = 0.3, (d) 5 — 0.4. The parameters of the unstable 
focus were set at A = 0.1 and w = tt as previously. It is 
observed that increasing the modulation amplitude e for 



a constant latency value generally leads to an enlarge- 
ment of the K interval of successful control, depending 
on the type of the modulation. For a sufficiently large 
latency time when the Pyragas control fails (e.g, 5 — 0.4 
in panel (d) of Figs. [T0llT2)) . inclusion of a variable delay 
in the controller makes the control possible again for a 
suitable choice of the modulation amplitude e > 0. Nev- 
ertheless, the control domain depends on the modulation 
type, and it is seen that in this case it is much less pro- 
nounced for a sawtooth-wave modulation than for the 
other two delay-modulation types. 



In a similar fashion, from Eqs. p7p and (|55)) one 
obtains the parametric representation of the stability 



10 




- u 


i. 


--0.3 


0.75 


-0.6 


S 0.5 


-0.9 


0.25 


-1.2 

n ^' 


0. 
1. 


- -0.3 


0.75 


-0.6 


S 0.5 


-0.9 


0.25 


-1.2 


0. 



FIG. 9: (Color online) Same as Fig. [8] for a sine- wave modu- 
lation of the delay for fixed e — 0.5 and different values of the 
latency time: (a) 5 = 0.2, (b) S — 0.4, and for a square-wave 
modulation: (c) S — 0.2, (d) 5 = 0.4. Other parameters are 
the same as in Fig. [8] 




0.75 


\ 


(c) 


S 0.5 


/^ 


- 


0.25 

A 


/ 


- 



1 



6 ( 




K 



FIG. 10: (Color online) Control domain in the {K, e) plane for 
a sawtooth-wave modulation of the delay for different values 
of the latency time: (a) S = 0.1, (b) S = 0.2, (c) 5 — 0.3 and 
(d) 5 = 0.4. The nominal time delay is fixed at to — 1. The 
parameters of the free-running system are as in Fig. (2] 



boundaries in the {K, e) plane 

A sin n{S + tq) + {u — il) cos n{d + tq) 






sin(J7ro) 



A sin(r2(5) + {io -fl) cos{QS) 



(43) 



A sin 17(5 + To) + {uj — ^)cosQ{6 + tq) ' 



(44) 



in which the dependence of e on the cigenfrequency fJ 
enters implicitly in the second equation via the function 
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FIG. 11: (Color online) Control domain in the {K,e) plane 
for a sine- wave modulation of the delay for different values of 
the latency time: (a) 5 = 0.1, (b) S = 0.2, (c) 5 = 0.3 and (d) 
S = 0.4. The other parameters are as in Fig. 1101 




FIG. 12: (Color online) Control domain in the {K, e) plane 
for a square- wave modulation of the delay for different values 
of the latency time: (a) S = 0.1, (b) S = 0.2, (c) S — 0.3 and 
(d) S = 0.4. The other parameters are as in Fig. 1101 



X- The calculated boundaries of the control domain are 
the solid lines depicted in Figs. [TUHT^ 

Overall, we see a massive destruction of the control 
domains with increasing latency time. Compared to the 
effect of a phase rotation in the previous section, the 
presence of control-loop latency tends to spoil the con- 
trol mechanism radically. Besides the destructive inter- 
ference effect, on which variable-delay feedback control 
mainly relies, an instantaneous feedback provides the 
main source of stability which normally is reflected in 
an extensive coverage of the parameter space with so- 
lutions of successful control. If this dissipation term is 
replaced by a much less effective term due to latency, the 
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control performance is consequently lost. Thus in any 
experimental application of variable-delay feedback con- 
trol one should split the control term if possible, so that 
the instantaneous part is implemented by a direct mod- 
ification of the system to be controlled, and the variable 
delay part can be realized by separate devices. 



V. CONCLUSIONS 

Extension of the Pyragas time-delayed feedback con- 
trol of unstable equilibria by introducing time- varying de- 
lay in the delayed term of the feedback force significantly 
enlarges the stability domain in the control parameter 
space. 

We have shown that an analytical investigation of the 
control domains becomes possible in the range of high- 
frequency delay modulations, in which case the variable- 
delay term in the controller can be approximated by an 
equivalent distributed delay term. Restricting our atten- 
tion to this high-frequency limit, we have investigated the 
control efficiency of the variable-delay feedback method 
by considering a simple normal form of an unstable fo- 
cus and two different realizations of the controller which 
are experimentally relevant: a non-diagonal coupling of 



the control force realized via a phase-dependent coupling 
matrix, and the incorporation of an additional small con- 
stant delay term in all the arguments of the feedback 
force that represents the latency of the control-line real- 
ization. In each case, the variable-delay feedback control 
with a finite modulation of the delay is shown superior 
with respect to the constant delay case of Pyragas con- 
trol. We have also shown that, as in the original time- 
delayed feedback control, the torsion of the orbit is a 
necessary condition for the control method to be able to 
stabilize equilibria. 

Since in practice the equivalence of the distributed de- 
lay case holds already for fairly low modulation frequen- 
cies, we propose that variable-delay feedback is a conve- 
nient experimental method for realizing distributed delay 
feedback with different delay distribution. 
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